Confinement of Skyrmion states in noncentrosymmetric magnets 
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Skyrmionic states in noncentrosymmetric magnets with Lifshitz invariants are investigated within 
the phenomenological Dzyaloshinskii model. The interaction between the chiral Skyrmions, being 
repulsive in a broad temperature range, changes into attraction at high temperatures. This leads 
to a remarkable confinement effect: near the ordering temperature Skyrmions exist only as bound 
states, and Skyrmion lattices are formed by an unusual instability-type nucleation transition. Nu- 
merical investigations on two-dimensional models demonstrate the confinement and the occurrence 
of different Skyrmion lattice precursor states near the ordering transition that can become thermo- 
dynamically stable by anisotropy or longitudinal softness in cubic helimagnets. We introduce a new 
fundamental parameter, confinement temperature T p separating the peculiar region with "confined" 
chiral modulations from the main part of the phase diagram with regular helical and Skyrmion 
states. 
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In certain noncentrosymmetric magnetic systems, the 
asymmetric Dzyaloshinskii-Moriya (DM) exchange re- 
sults in chiral couplings that can stabilize long-period, 
non-collinear modulations of the magnetization with a 
fixed sense of rotation [TJ [2]. These chiral couplings 
are phenomenologically described by Lifshitz invariants 
that destroy the homogeneity of ordered phases pQ. In 
highly symmetric systems with Lifshitz invariants multi- 
ple modulations occur as textures with localized twists in 
two or more spatial directions |3J Hj • In noncentrosym- 
metric magnetic systems with Lifshitz invariants these 
double-twist motifs exist as smooth (baby)-Skyrmions 
[5HZ], static solitonic textures localized in two spatial 
directions, which can be extended into the third direc- 
tion as Skyrmion strings or Hopfions. These magnetic 
Skyrmions are stabilized solely by the chiral DM cou- 
plings [7] , which prevent a spontaneous collapse into 
topological singularities. Skyrmionic matter is created 
by the condensation of these solitonic units, similar to 
vortex matter in type-II superconductors [5]. Just such 
chiral Skyrmions have been recently observed in thin 
layers of noncentrosymmetric ferromagnet (Fe,Co)Si [5]. 
Skyrmionic states stabilized by Lifshitz-type invariants 
may exist and form extended mesophases in various con- 
densed matter systems, as chiral liquid crystals, ferro- 
electrics, multiferroics, and in confined achiral systems 
(e.g., thin magnetic layers) [HI5HTT]. 

Skyrmionic textures usually form through nucleation, 
following a classification introduced by DeGennes [T2] 
for (continuous) transitions into incommensurate mod- 
ulated phases. As demonstrated in Refs. 0[6], for the 
low-temperature micromagnetic model of chiral magnets, 
at the transition from the field-driven Skyrmion lattices 
into the polarized homogeneous state the lattice period 
diverges and the Skyrmions are set free as localized exci- 
tations. As the Skyrmions retain their size and axisym- 
metric shape, there is a full spectrum of lattice modes up 
to the transition, in contrary to an instability type tran- 



sition where the amplitude of one fundamental mode acts 
as small parameter of the transition [12 . 

Here, we show for the standard model of chiral 
isotropic ferromagnets [TJ [7J [T3] that Skyrmions are con- 
fined very close to the ordering temperature. In that 
part of the phase diagram, the creation of Skyrmions as 
stable units and their condensation into extended tex- 
tures occurs simultaneously through a rare case of an 
instability-type nucleation transition |14j . but the con- 
fined Skyrmions as discernible units may arrange in dif- 
ferent mesophases. This is a consequence of the coupling 
between the magnitude and the angular part of the order 
parameter. Thus, near the ordering transitions, the local 
magnetization is not only multiply twisted but also lon- 
gitudinally modulated. From numerical investigations on 
2D models of isotropic chiral ferromagnets, a staggered 
half-Skyrmion square lattice at zero and low fields and a 
hexagonal Skyrmion lattice at larger fields are found in 
overlapping regions of the phase diagram near the tran- 
sition temperature. Furthermore, the thermodynamic 
stability of Skyrmionic states can be favoured with re- 
spect to one-dimensional modulations by supplementing 
the model with cubic exchange anisotropy or in a modi- 
fied model for metallic chiral magnets [7]. 

Within the standard phenomenological (Dzyaloshin- 
skii) theory [T] the magnetic energy density of a non- 
centrosymmetric ferromagnet can be written in the di- 
mensionless form [Jj [13] 

<f> = (grad m) 2 — iod(ixi) — h(n ■ m) + am 2 + m A . (1) 

Here, reduced values of the spatial variable x = r/Ljj, 
the magnetization m = M/Mq, and the applied magnetic 
field hn (h = \H\/Ho, n is a unity vector along H), are 
expressed via the parameters of the energy density 1 13] 
w(M) = A(gradM) 2 + Dw D (M) - H.-M + f(M) (where 
f(M) = aiM 2 + a 2 Af 4 , and M = |M|): L D = A/D, 
H = kM , M = (K/a 2 ) 1/2 , a = ai/«, k = D 2 /(A). 




FIG. 1: (Color online) (a) Equilibrium solutions of local- 
ized Skyrmions, shown as cross section in the inset, display 
strong localization of 6(p) (blue) and weak variation of the 
magnetization modulus m(p) (red) in a broad range of tem- 
perature and field (a,h). Exemplary profiles for a = —0.25, 
h — 0.25. (b) At lower fields, isolated Skyrmions condense 
into a dense-packed hexagonal lattice, (c) Near the ordering 
temperature a c = 0.25 square lattice solutions with staggered 
half-Skyrmion cells arise. 

DM energy w£> consists of Lifshitz invariants [T] 

pfy = mi(drrij/dxk) - rrij(dmi/dxk). (2) 

For noncentrosymmetric uniaxial ferromagnets DM func- 
tional wd are listed in [5]. In particular, for isotropic 
and cubic helimagnets u>£>(m) = C^J + C z X y + Cxy 
— mrotm [13] (for details see appendix). Functional ([!]) 
includes only the basic (isotropic) interactions essential 
to stabilize chiral modulations. This depends on three 
internal variables (components of the magnetization vec- 
tor m) and two control parameters, the reduced magnetic 
field amplutude h and the "effective" temperature a(T). 
In chiral magnetism models ([I]) play a fundamental role 
and are similar to the Frank energy in liquid crystals and 
Landau- Ginzburg functional in superconductivity. We 
consider 2D chiral modulations homogeneous along the 
applied field h||z and modulated in the plane perpendic- 
ular h. 

Isolated and embedded Skyrmions. The equations min- 
imizing functional include solutions for axisymmet- 
ric localized states (isolated Skyrmions), ip — ip(4>), and 
0(p), m(p) where we use spherical coordinates for the 
magnetization (m, 9, yj) and cylindrical coordinates for 
the spatial variables (p,4>,z). The solutions "0(0) f° r all 
noncentrosymmetric classes have been derived in [5] . For 
cubic helimagnets and uniaxial systems with n22 sym- 
metry, i/j — (jy — 7r/2 (Fig. [l]a). Profiles 9(p) and m(p) of 
isolated Skyrmions are derived from numerical solutions 
of the Euler equations, as given in Ref. Hand common for 
all noncentrosymmetric classes with Lifshitz invariants, 
with the boundary conditions 8(0) = n, dm/dp(0) = 0, 
9(oo) = 0, to(oo) = mo (mo is the magnetization in 
the saturated state). For extended textures of two- 
dimensional models, the functional has been investigated 
by numerical energy minimization using finite-difference 
discretization on rectangular grids with adjustable grid 
spacings [7]. 

Isolated Skyrmions (Fig. [I] (a)) exist only below a crit- 
ical line ho and condense into a hexagonal lattice (Fig. [2]) 
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FIG. 2: (Color online) Field h vs. temperature a phase 
diagram, (a) Areas with repulsive (I) (shaded), attractive 
(II) Skyrmions (hatched), and strictly confined pocket (III) 
(cross-hatched) are separated by h* (dashed line), Eq. J3J. 
Isolated Skyrmions collapse at critical line ho- Above this 
line no static Skyrmions exist. Below line h c Skyrmions con- 
dense into a hexagonal lattice In region (II), the hexagonal 
Skyrmion lattice (SL) exists as metastable state up to the 
nucleation field h n . For temperatures between points A and 
B, h n < ho, for larger temperatures as < a < a c — 0.25 the 
isolated Skyrmions disappear at lower fields than the dense 
SL, ho > h n . For clarity, line h n is only schematically given in 
panel (a), numerically exact data are shown in panel (b). De- 
tail of the phase diagram (c) near the ordering temperature 
shows the existence regions for different modulated states. 
Lines for first order transitions: a square half-SL o hexago- 
nal SL, p helicoid <-» hexagonal SL. Line 7 marks the continu- 
ous transition from the conical equilibrium phase in isotropic 
systems to the paramagnetic phase. 

below a field h c , which marks the transition between the 
Skyrmion lattice (SL) and the homogeneous paramag- 
netic state (see Fig. [2] (a), (b)). Near the ordering tem- 
perature a square half-Skyrmion lattice (Fig. [I] (c)) has 
lower energy than the hexagonal lattice. Half-SLs con- 
sist of cells with up and down magnetization in the cen- 
ter and in-plane magnetization along the cell boundaries. 
Such cells have a topological charge 1/2. In the hexag- 
onal SLs the magnetization at the boundaries (center) 
is (anti)parallel to the applied field. The cells bear unit 
topological charges. 

Confinement. By solving the linearized Euler equa- 
tions for the asymptotics of isolated Skyrmions Am = 
(m — TOq), cx exp(— up) (p 3> 1) one finds three distinct 
regions in the magnetic phase diagram with different 
character of Skyrmion- Skyrmion interactions (Fig. [2] (a)): 
repulsive interactions in a broad temperature range (area 
(I)) are changed to attractive interaction at higher tem- 
peratures (area (II)). Finally in area (III) near the order- 
ing temperature a c — 0.25 strictly confined Skyrmions 
exist. These regions are separated by the line 

h* = y/2±P(a)(a + 1 ± P(a)/2), P(a) = V3 + 4a (3) 

with turning points p (-0.75,-^/2/4), q (0.06,0.032^/5), 
and u (-0.5, 0) (dashed line in Fig. [2] (a)). In the major 
part of the phase diagram, the evolution of SLs under a 
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FIG. 3: (Color online) Numerically exact solutions for hexag- 
onal (a)-(c) and square (d)-(f) Skyrmion lattices within the 
confinement pocket. Contour plots of M z (x,y) and M(x,y) 
(a),(d) and their diagonal cross-sections (b),( e), (c),(f) are 
derived by minimization of $ with a — 0.23 and different 
values of the applied field. Solid (dashed) lines in panels (e), 
(f) indicate profiles for cells with (anti)parallel magnetization 
in the center. 



magnetic field closely agrees with the behavior studied 
earlier for the low-temperature limit jSJ [5] and the tran- 
sition mechanism at the high-field limit is of the nucle- 
ation type with isolated Skyrmion excitations appearing 
below the instability line ho. The temperature of the 
turning point p (a p = —0.75), a confinement tempera- 
ture T c determines the frontier between low- (T < T p ) 
and high-temperature (T p < a < T c ) chiral modulations. 
Usually AT = T c — T p CT C and high-temperature mod- 
ulations are restricted by a narrow vicinity of the Neel 
temperature. According to calculations in [7] for MnSi 
high-temperature Skyrmions exist in temperature inter- 
val AT = T c — T p — 2 K below the Curie temperature 
T c . 

Evolution of confined Skyrmions. In the vicinity of the 
ordering temperature the Skyrmion lattices exist within 
the confinement pocket (III). In this region Skyrmion 
states drastically differ from those in the main part of 
the phase diagram. Due to the "softness" of the magne- 



tization modulus the field-driven transformation of the 
Skyrmion lattices evolves by distortions of the modulus 
profiles m(p) both in the hexagonal and square Skyrmion 
lattices while the equilibrium periods of the lattices do 
not change strongly with increasing applied field (Fig. 
[3]). Despite the strong transformation of their internal 
structures the Skyrmion lattices preserve axisymmetric 
distribution of the magnetization near the centers of the 
Skyrmion lattice cells (Fig. [3|. This remarkable property 
reflects the basic physical mechanism underlying the for- 
mation of Skyrmion lattices. The local energetic advan- 
tage of Skyrmion lattices over helicoids is due to a larger 
energy reduction in the "double-twisted" Skyrmion cell 
core compared to "single- twisted" helical states [4j [7]. 
This explains the unusual axial symmetry of the cell cores 
and their stability. An increasing magnetic field gradu- 
ally suppresses the antiparallel magnetization in the cell 
core reducing the energetic advantage of the "double- 
twist" and increases the overall energy of the condensed 
Skyrmion lattice. At critical field h c a first-order tran- 
sition occurs into the polarized paramagnetic state. In 
the interval h c < h < h n the hexagonal Skyrmion lattice 
exists as a metastable state. At the lability field h n the 
magnetization modulus in the cell center becomes zero 
(see magnetization profile for h = 0.042 in Fig. [3] (c)). 
At this field, the " double- twist" region is suppressed, and 
the lattice loses its stability. 

At zero field with increasing temperature a the mag- 
netization modulus m in hexagonal and square lattices 
gradually decreases to zero at the ordering temperature. 
This is the instability-type nucleation transition into the 
paramagnetic phase : the order parameter m becomes 
zero at the transition point (as in the instability mode), 
however, the lattices retain their symmetry and the ar- 
rangement of axisymmetric Skyrmions up to the criti- 
cal point. In finite fields in the region (II) the attrac- 
tive Skyrmion-Skyrmion interaction means that multi- 
Skyrmions can always form bound states. Near a c such 
clusters of Skyrmions are more stable than the isolated 
Skyrmions. This is seen in Fig. [2] (a),(b). For tempera- 
tures above point B, as < a < a c , the metastable SL as 
the densest infinite cluster is more stable than isolated 
Skyrmions. Thus, in region (II) of the phase diagram, 
there is clustering and, for higher temperatures, confine- 
ment of isolated Skyrmion excitations. 

The predicted Skyrmionic textures and the pre- 
cursor phenomena associated with the confinement of 
Skyrmions near the ordering transition are observable 
in magnetic materials with appropriate symmetry. In a 
large group of uniaxial noncentro symmetric magnets the 
DM energy is described by gradients only along directions 
perpendicular to the axis (e.g. multiferroic BiFeOs, space 
group i?3c and antiferromagnetic Ba2CuGe207, space 
group PA2 1 m [HE3). In such magnets one-dimensional 
modulations with the propagation vector in the basal 
plane [helicoids) exist in broad ranges of the magnetic 



4 




0.21 0.23 0.25 0.18 0.22 0.26 



FIG. 4: (Color online) Magnetic phase diagrams of cubic heli- 
magnets with exchange anisotropy b = —0.05 and the applied 
field along (111) (solid) and (001) (dashed) axes (a) contains 
regions with thermodynamically stable hexagonal SLs and 
half-SLs. Magnetic phase diagram of the modified isotropic 
model with n — 0.8 (b) has extended domains with thermo- 
dynamically stable helicoid, half-SLs and hexagonal SLs with 
the magnetization along the cell axes parallel to the magnetic 
field (hatched area). 

fields [3 [T5] . In the confinement pocket of the phase di- 
agram (III) the helicoids also exist only as bound states 
[US] , They have the lowest energy at lower fields, and the 
hexagonal SL corresponds to the global energy minimum 
at higher fields (Fig. [2] (c)). 

Cubic helimagnets. In other groups of chiral mag- 
nets the DM energy includes contributions with gradients 
along all three spatial directions. This stabilizes chiral 
modulations with propagation along the direction of an 
applied field as cone phases [T3] . For the isotropic model 
<i>(m) ([!]) the cone phase solution with the fixed magne- 
tization modulus and rotation of m around the applied 
magnetic field: if> = z, cos(#) = h/m, m = \a — l/2|/2, 
is the global energy minimum in the whole region where 
modulated states exist. For cubic helimagnets, the en- 
ergy density ([I]) has to be supplemented by anisotropic 
contributions, $ a = bj^iidnii/dxi) 2 + kJ2i m i7 where 
b and k are reduced values of exchange and cubic 
anisotropies |13| . These anisotropic interactions impair 
the ideal harmonic twisting of the cone phase and lead 
to the thermodynamic stability of Skyrmion states, as 
shown in the equilibrium phase diagram Fig. |i](a). The 
difference between the energy of the hexagonal Skyrmion 
lattice W s k and of the cone phase W cone calculated for 
the isotropic model, AW m i n = W s k — W cone , has min- 
ima along a curve h x (a) which reaches the critical point 
h K (a c ) = as AWmin = 0.0784(0.25 - a)). Weak ex- 
change anisotropy of a cubic helimagnet, therefore, cre- 
ates a pocket around a c , where the hexagonal Skyrmion 
lattice becomes the global energy minimum in a field 
(Fig. [1] (a)). This case is realized in cubic helimag- 
nets with negative exchange anisotropy (6 < 0) as in 
MnSi [13] . This anisotropy effect provides a basic mech- 
anism, by which a Skyrmionic texture is stabilized in 
applied fields, as observed experimentally as so-called 
"A-phase" in MnSi [T7HH]- The exchange anisotropy 



b < also leads to the thermodynamic stability of half- 
SLs (Fig. [4] (a)). The stabilization of such textures may 
be responsible for anomalous precursor effects in cubic 
helimagnets in zero field [7J [501 121] • 

Chiral modulations in non-Heisenberg models. A gen- 
eralization of isotropic chiral magnets proposed in [7J 
replaces the usual Heisenberg-like exchange model by 
a non-linear sigma-model coupled to a modulus field 
with different stiffnesses. This yields a generalized gra- 
dient energy for a chiral isotropic system with a vec- 
tor order parameter, which is equivalent to the phe- 
nomenological theory in the director formalism [4j [Jj 

Eijto"*) 2 -> Ei,#i0 2 + (1 - ^EijCft"*) 3 = 

m 2 J2i j(9i n j) 2 + vJ2i(di m ) 2 ■ Parameter n equals unity 
for a "Heisenberg" model, in chiral nematics 77 = 1/3 0J. 
Confined chiral modulations are very sensitive to values 
of rj < 1. The magnetic phase diagram calculated for 
rj = 0.8 includes pockets with square half-SLs, hexagonal 
SLs with the magnetization in the center of the cells par- 
allel to the applied field, and helicoids with propagation 
transverse to the applied field (Fig. |I]b). 

The confinement effects on chiral Skyrmions strongly 
changes the picture of the formation and evolution of chi- 
ral modulated textures and shed new light on the prob- 
lem of precursor states observed as blue phases in chiral 
nematics [3] and in chiral magnets [7J H7T121) . The re- 
sults show that confinement / dcconfincmcnt transitions 
of localized string-like solitons can be realized in these 
condensed matter systems. They provide counterparts 
of formation mechanisms for extended microscopic mat- 
ter from topological solitons, as devised originally in the 
Skyrme model [22 . 
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APPENDIX:DZYALOSHINSKII THEORY OF 
CHIRAL HELIMAGNETISM 

The appendix includes technical comments on the 
mathematical structure of the interaction functional |l]) , 
on the mathematical relation between this model and 
those for different classes of noncentrosymmetric mag- 
netic systems and chiral liquid crystals. The equations 
for isolated and bound Skyrmions are provided along 
with a short information on the numerical methods to 
solve these problems. 

1. Basic model. 

In noncentrosymmetric ferromagnets chiral asymmetry 
of exchange coupling related to the quantum-mechanical 
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Dzyaloshinskii-Moriya interactions [23] induces long- 
period modulations of the magnetization vector M with 
a fixed rotation sense [1] . The equilibrium chiral modula- 
tions are derived by minimization of the magnetic energy 
which can be written in the following general form pQ 

W = A^2(dM ] /dx l ) 2 -M-H + W (M) + W D {M) (4) 

where A is the exchange stiffness constant, H is an ap- 
plied magnetic field, Wo collects short-range magnetic in- 
teractions independent on spatial derivatives of the mag- 
netization, and Dzyaloshinsky-Moriya energy Wd is com- 
posed of antisymmetric terms linear with respect to first 
spatial derivatives of the magnetization vector M(r) 

4f = M t {BMj/dxk) - M % (dMj/dxk) . 

Functional forms of Wd energy contributions are deter- 
mined by crystallographic symmetry of a noncentrosym- 
metric magnetic crystal [T] [5] [T3] . Particularly, for im- 
portant cases of noncentrosymmetric ferromagnets be- 
longing to cubic (23, 432) and uniaxial crystallographic 
classes (nmm, 42m) the DM energy contributions have 
the following form [U [5j [13] 

W D = D OCW + 4*) + £<fj) = DM • rotM 

for (23), (432); 
Wd — D (£ W + £(») ) for (nmm) ; 

W D =D(£<g+C$) for (42m) (5) 

where n = 3,4,6 and D is a Dzyaloshinskii constant. 
For other noncentrosymmetric classes this Wd energy in- 
cludes two or more Dzyaloshinskii constants [5]. 

Near the ordering temperatures the magnetization am- 
plitude varies under the influence of the applied field and 
temperature. Commonly this process is described by in- 
cluding into the magnetic energy an additional term [13] 

W -> W + a x M 2 + a 2 M 4 , a x = J(T - T c ), a 2 > (6) 

This model is fundamental to the systematic phe- 
nomenological description of magnetic states in noncen- 
trosymmetric magnetic systems pQ. Interaction func- 
tional Q plays in chiral magnetism a similar role as 
the Frank energy in liquid crystals [24] or Ginzburg- 
Landau functional in physics of superconductivity |25j . 
Its form is dictated by the natural extension of Landau's 
approach to ordering transtions for modulated systems 
which break the third Lifshitz condition, as pioneered by 
Dzyaloshinskii in particular in magnetism [26j . As a re- 
sult, this type of theory constitutes the canonical form 
of statistical field theories for condensed matter systems, 
where particular couplings and symmetry enable Lifshitz 
invariants to induce modulated phases. Extensions in- 
cluding higher-order secondary effects in magnetism (like 
anisotropies or dipolar interactions) are widely used to 



describe evolution of modulated states in many chiral 
magnets [6] [T] [27] , including metamagnetic transitions 
in cubic helimagnets induced by high pressure [T3 ] l2"8 l l2"9 ] . 

By rescaling the spatial variable in |l]),((6| x = r/L D , 
the magnetic field h = H/ifo, the magnetization m = 
M/M 

L D = A/D, H = kM , Mo = («/a 2 ) 1/2 , 

K = D 2 /(2A), a = oi/k = J(T - T c )/k. (7) 

energy W f6j) can be written in the following reduced 
form (Eq. Q) 

$ = (grad m) 2 — w/)(m) — h(n ■ m) + am 2 + m 4 , (8) 

where h = |h| and n is a unity vector along the ap- 
plied magnetic field. Functional ^ includes three inter- 
nal variables (components of the magnetization vector 
m) and two control parameters, the reduced magnetic 
field amplutude h and the "effective" temperature a(T) 
([7]). By direct minimization of Eq. ^ we have derived 
one-dimensional (kinks, helicoids, conical helices or spi- 
rals) and two-dimensional solutions (isolated and bound 
Skyrmions) for arbitrary values of the control parame- 
ters. These results are collected in the phase diagram of 
solutions (Fig. 2). 

Energy (JsJ includes only basic interactions essensial 
to stabilize Skyrmionic and helical phases. Solutions 
for chiral modulated phases and their most general fea- 
tures attributed to all chiral ferromagnets are deter- 
mined by interactions functional dsh. Generically, there 
are only small energy differences between these differ- 
ent modulated states. On the other hand, weaker en- 
ergy contributions (as magnetic anisotropy, stray-fields, 
magneto-elastic couplings) impose distortions on solu- 
tions of model ^ which reflect crystallographic symme- 
try and values of magnetic interactions in individual chi- 
ral magnets. These weaker interactions determine the 
stability limits of the different modulated states (Fig. 
4). The fact that thermodynamical stability of individual 
phases and conditions of phase transfomations between 
them are determined by magnetocrystalline anisotropy 
and other relativistic or weaker interactions means that 
(i) the basic theory only determines a set of different 
and unusual modulated phases, while (ii) the transitions 
between these modulated states, and their evolution in 
magnetization processes depends on symmetry and de- 
tails of magnetic secondary effects in chiral magnets, in 
particular the strengths of relativistic magnetic interac- 
tions. Thus functional ^ is the generic model for a man- 
ifold of interaction functionals describing different groups 
of noncentrosymmetric magnetic crystals, because it al- 
lows to identify the basic modulated structures that may 
be found in them. 
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2. Equations for isolated and bound Skyrmions 

The Euler equations for functional Q have solutions 
for axisymmetric isolated structures of type tp = ip((j>), 
0(p), m(p) (we introduce here spherical coordinates for 
the magnetization (m, 9, ip) and cylindrical coordinates 
for the spatial variables (p, <p, z)) 7J. Solutions ip = ■0(0) 
are known for all uniaxial and cubic noncentrosymmet- 
ric ferromagnets [5]. Particularly, ip = <p — tt/2 for cu- 
bic helimagnets, ip = <p and ip = — (p — n/2 for uniax- 
ial ferromagnets with (nmm) and 42m symmetry, corre- 
spondingly. After substitution of solutions ip — ip((f>) 
into ([Sjl and integration with respect to cf> the total 
energy E of a Skyrmion (per unit length along z) is 
E = 2tt J °° $(m, 8)pdp with energy density 
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where a common convention df/dx = f x is applied. The 
Euler equations for the functional ^ 
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with boundary conditions 9(0) = it, 9(oo) = 0, m(oo) = 
mo describe the structure of isolated Skyrmions (the 
magnetization of the homogeneous phase mo is derived 
from equation 2am + 4m 3 — h — 0). Typical solutions 
9(p), m(p) of Eqs. (10) are plotted in Fig. 1. 



For present work, we have evaluated numerically the 
solutions for the energy functionals of type Eq. within 
a standard approach using finite differences for gradi- 
ent terms and adjustable grids to accommmodate mod- 
ulated states with periodic boundary conditions. Search 
for modulated states and energy minimization was done 
by Monte Carlo simulated annealing. For 2-dimensional 
models, this approach can give converged results, as 
checked by comparison with analytical results. Including 
secondary effects, the approach is able and has been used 
by us to calculate thermodynamic stability of phases, 
and to study the complete temperature-field phase di- 
agrams in terms of equilbrium (mean-field) states of the 
energy functionals Eq. ([!]). As we are investigating here 



long-period modulated phases, the phase diagrams are 
expected to give qualitatively the correct results for these 
phenomenological models. 
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